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ABSTRACT 

Helioseismic investigation has suggested to apply the turbulent convection models (TCMs) to the convective 
overshooting. Using the turbulent velocity in the overshooting region determined by the TCM, one can deal 
with the overshooting mixing as a diffusion process, which leads to incomplete mixing. It has been found that 
this treatment can improve the solar sound speed and the Li depletion in open clusters. In order to investigate 
whether the TCM can be applied to the overshooting mixing outside the stellar convective core, new observa- 
tions of the eclipsing binary star HY Vir are adopted to calibrate the overshooting mixing parameter. The main 
conclusions are as follows: (i) the TCM parameters and the overshooting mixing parameter in the solar case are 
also suitable for the eclipsing binary system HY Vir; (ii) the incomplete mixing results in a continuous profile 
of the hydrogen abundance; (iii) the e-folding length of the region, in which the hydrogen abundance changes 
due to the overshooting mixing, increases during the stellar evolution. 
Subject headings: convection — turbulence — stars: individual: HY Vir 



1. INTRODUCTION 

The convective overshooting in the stellar interior is caused 
by turbulent flows moving across the standard convective 
boundary, which is defined by the Schwarzchild criterion (i.e., 
Vji > V a d), into the dynamically stable region. This pro- 
cess leads to chemical mixing in the overshooting region. 
Based on the framework of non-local mixing length theo- 
ries, the overshooting region is considered as a fully mixed 
region with a length Iqv = olovHp, where Hp is the lo- 
cal pressure scale height and aov is a parameter. But the 
helioseismic investigation has pointed out that only turbulent 
convection model s (TCMs) can properly describe the convec- 
tive overshooting (Christe nsen-Dalsgaard et aTJl20 1 11) . There- 
fore, it is a tendency to apply TCMs to the convective over- 
shooting in the stellar interior. TCMs are non-local convec- 
tion theories based o n the hydrodynamic equa t ions and mod- 
eling assumptions dXiond|198lL fl985L 11989b ICanutol 1 1 9971; 
Canuto & Dubo vikovll 19981; IXiong & Deng||200U iLi & Yang) 
200 U iDeng et all 120061; ILi & Yandl2007tTCanutoll201 U iLil 
2012). Using the turbulent velocity determined by TCMs, one 
can deal with the overshooting mixing as a diffusion process 
which results in incompl ete mixing in the overs hooting re- 
gion IZhang & Lil (I2012al) have applied the TCM (ILi & Yangj 
2007) and the incomplete overshooting mixing to the down- 
ward overshooting region of the solar convective envelope and 
found that the sound speed profile is significant improved. 
Applying the same method to the stellar model also shows 
improvements on the Li depletion in open cluster, thus the 
TCM has been considered to be generally applicable to the 
downward over shooting regio n of the convective envelope for 
low-mass stars (Zhang 2 012b . 

It is interesting whether the incomplete overshooting mix- 
ing based on the TCM is suitable for the convective core over- 
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shooting. The convective core overshooting refreshes hydro- 
gen in the nuclear burning core, deeply affecting the evolu- 
tions of massive and intermediate-mass stars. The evolution 
track of a star with a convective core in the Hertzsprung- 
Russel diagram is sensitive to the efficiency of the core over- 
shooting mixing. Accordingly, it is a direct and efficient way 
of studying the core overshooting mixing by calibrating the 
radius and effective temperature of the stellar model to fit 
the observations. Double lined detached eclipsing binary sys- 
tems in the main sequence phase are the best candidates to 
do that, because the mass, radius, and temperature of each 
comp onent of a b i nary s yste m can be obta ined by observa- 
tions. iRibas et all (120001) and lClaretl (120071) have studied the 
core overshooting based on a set of double lined eclipsing 
binary systems and found that the fully mixed overshooting 
region with a length 0.20 ~ 0.25iJp is the overall best. Re- 
cently, Lacv & Fekelj d201 lb have observed the double lined 
detached eclipsing binary system HY Vir and acquired very 
accurate masses and radii of the two components from the 
analysis of new light curves and radial velocity curves. Th e 
mass of each component of HY Vir is ( Lacy & Feke |i20lTl) : 
Mi = 1.838 ± O.OO9M and M 2 = 1.404 ± 0.006M Q for 
the primary and the secondary, respectively. These results in- 
dicate that both the primary and the secondary have the con- 
vective core. The accurate data make HY Vir a good candidate 
to study the core overshooting mixing described by the TCM. 

In this paper, I investigate the suitable parameter of the 
overshooting mixing via calibrating the effective temperatures 
and the radii of the stellar models of HY Vir to fit the obser- 
vations, and stud y the effects of the application of the TCM 
dLi & Yang||2007b . The method of applying the TCM in the 
stellar models is described in Section 2. The calibrated results 
for the system of HY Vir are shown in Section 3. The proper- 
ties of the stellar models based on the TCM are discussed in 
Section 4. The conclusions and discussions are presented in 
Section 5. 



2. THE METHOD OF CALCULATION 
2.1. The overshooting mixing 
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In the traditional way, the overshooting is assumed to re- 
sult in the completely mixed region with the length being 
proportional to the local pressure scale height. This de- 
scription of the overshooting is based on the framework of 
the 'ballistic' overshoot models, which are non-l o cal mix- 
ing length theories(e.g . . iShayiv & Salpepteri (11973 ); Maeder 
d!975l) : iBressan et all dl981l) ; IZahnl dl991l) ). Those non- 
local mixing length theories result in a steep temperature 
gradient p rofile at the boundary of the oversh ooting re- 
gion (e.g., iShaviv & Salpepteri (11973b ; IZahnl (119911) ). How- 
ever, the recent helioseismic investigation has found that 
the required temperature gradient in the overshooting re- 
gion is not in agreement with the framework of the 'bal- 
listic' overshoot models, but suggests that only TCMs are 
capabl e of providing required tem perature gradient profile 
dChristensen-Dalsgaard et aT]|201 11) . 

According to the hydrodynamic equations, the chemical 
mixing is determined by the turbulent velocity - turbulen t con- 
centration correlation dXiongll 1 9861; ICanutoll 1999L 120 1 ll) . Un- 
fortunately, there is a lack of the knowledge for such a correla- 
tion in the astrophysical case. Besides the fully mixing, the in- 
complete mixing described by the diffusion is at p resent pop- 
ular on de aling with the oversh o oting mixing (e.g..|Deng et alJ 
( 1996a lbh;lFrevtag et alJ(fT996h ; lVentura euSdl998h ; lHerw!g 
(2000); iPaxton et alJ d201 lb ). In the diffusion description, 
the turbulent velocity in the overshooting region is required 
to calculate the diffusion coefficient. The turbulent velocity 
is usually set as an expo nentially decreasing f unction in the 
oversh ooting region (e.g., Frevtag et al. ( 1996T): [Ventura et all 
(1998)), and can be calculated by defining the exponential in- 
dex and the initial value (i.e., the turbulent velocity at the con- 
vective boundary). An option to determine the initial value 
is extrapolati o n from the mixing length th eory (MLT) (e.g., 
IFrevtag et all d 19961) ; iVentura et al] (119981) ). However, the 
MLT is a local theory, while the non-local effect dominates 
near the convective boundary. Another option is to obtain the 
turbulent velocity by solving the TCM, which is non-local tur- 
bulent convection theory. 

The diffusion coefficient is proportional to the characteris- 
tic length of mixing and the characteristic velocity: D = flv, 
where / is the characteristic length, v is the characteristic tur- 
bulent velocity, and / is a diffusion parameter, v is assumed 
as \fk, where k is the turbulent kinetic energy. The diffusion 
parameter multiplying the characteristic length is assumed to 
be proportional to Hp, i.e., fl = CxHp. Therefore, the dif- 
fusion coefficient of the overshooting mixing is assumed as 
follow: 



D ov = C x HpVh 



(1) 



The characteristic length of mixing I is between the Kol- 
mogorov scale and the largest eddy scale. In order to cause 
incomplete mixing in the overshooting region in a time scale 
comparable with the evolutionary time scale, the characteris- 
tic length o f mixing should be much smaller than the largest 
eddy s cale dDeng et al.|[l9 96a). By analogy with Eq.(27) & 
(29) in iDeng et al.1 d!996al) . CV in Eq.(l) should be on the 
magnitude order of 1 ~ 1 , which is suitable for the downward 
over shooting region of t he convective envelope in low-mass 
star dZhang & Lij|2012atlZhangll2012l) . 

The diffusion equation of the hydrogen abundance in the 
stellar interior is as follow: 



dt 



{{ATrpr 2 ) 2 D- 



where (OmixX / dt) means the time derivative of hydrogen 
abundance caused by only convective and overshooting mix- 
ing, and D is the diffusion coefficient due to convection or 
overshooting: 



D 



Dov;(Vr < Vad) 
Dcz] (Vfl > V ad ) 



(3) 



where Dcz is diffusion coefficient in the convection zone. 
In order to ensure fully mixing in the convection zone, large 
Dcz is required so that Dcz = Hpy/k is adopted in the 
calculations. It is not difficult to solve Eq.(2) by using the 
tridiagonal matrix algorithm. 

2.2. The stellar evolution code 

The s tellar evolu t ion c ode, which is originally devel- 
oped by iPaczynskil d 19691) and Kozlowski and updated by 
Sienkiewicz, is used to m odel the componen ts of HY Vir. The 
OPAL equation of state (Rog ers et al.ll 19961). the OPAL opac - 
ity tables for high temperatures (Iglesias & Rogers! 119961) . 
and the Alexander's opacit y tables for low temperatures 
dAlexander & Fergusol 119941) are used. The composition 
mixture is assumed to b e the same as the solar mixture 
dGrevesse & Sauvallfl998l) . 

In this paper, I solve the complete structure equations 
comprising the stellar structure equations and the TCM 
dLi & Y ang 2007) equations. In the stellar structure equations, 
the temperature gradient V is calculated as follow: 



V = Vb- 



T X Ur 



(4) 



dM r 



dM r ' 



(2) 



where u' r T' is the turbulent heat flux determined by the TCM. 

The evolution code has been modified to take into account 
the overshooting mixing. In the original code, the evolution 
of chemical composition in the stellar interior and the stellar 
structure are calculated separately. On each time step, the evo- 
lution of chemical composition (nuclear burning, complete 
mixing in the convection zone) is calculated before solving 
the stellar structure equations. In the present work, the proce- 
dure of solving Eq.(2) is performed at first. After the evolution 
of chemical composition (nuclear burning and the convective 
and overshooting mixing) during the time step being accom- 
plished, the complete structure equations are solved, and the 
stellar structure variables and the turbulent variables are up- 
dated. The scheme of solving the complete structure equa- 
tions is as follows: 

(1) Solve the TCM equations based on the current stellar 
structure variables (e.g., p, T, r, L). Calculate the temperature 
gradient V at all mesh points according to Eq.(4). 

(2) Solve the localized TCM in which the non-local terms 
(i.e., the diffusion terms) are ignored. Calculate the local 
temperature gradient Vi at all mesh points. Some discus- 
si ons about the local ized TCM can be found in Appendix A 
in lZhang&Lild2012bl) . 

(3) Calculate the ratio r\ = V / V l at all mesh points. Use r\ 
to get rf at all mesh points based on the relaxation technic, i.e., 
V' = Vpre + ^(V ~ Vpre), with a proper relaxation parameter 
£, where jf pre and rj pre are the previous value of rj and ?/. 

(4) Solve the stellar structure equations in which the tem- 
perature gradient is calculated as V — ij'Vi, and update stel- 
lar structure variables (e.g., p, T, r, L). It should be noticed 
in this process that Vl is updated when the stellar structure 
variables are updated in the Newtonian iterations. 
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TABLE 1 

The turbulent parameters adopted in this paper. 



c t 


C e 


C a 


Cti C e i Cfc a 


7.5 


0.2 


0.08 


2.5 0.85 



(5) Check differences of 1 77 — rf\ and \rj — 77 pre |. The cal- 
culations are thought to converge if \q — r/| and \r) — rj pre \ 
are less than required accuracy at all mesh points, otherwise 
return to step (1). 

It is not difficult to understand that the turbulent variables 
and the stellar structure variables are consistent with each 
other when T) pre = i) = if, thus the scheme is reasonable. 
The scheme is an i mproved version of the previous one (e.g., 
Zhang "&Ll (120091) 1. It is found that this scheme is more sta- 
ble and more efficient than the previous one; 

There are seven parameters in the TCM dLi & Ya ng 2007): 
a,Ct,C e are dissipation parameters of turbulent kinetic en- 
ergy k, turbulent heat flux u' r T', and turbulent temperature 
fluctuation T'T', respectively; C s , Cti, Cei are diffusion pa- 
rameters of k, u' r T' and T'T', respectively; Ck is the param- 
eter of 'return to isotropy' term. 

The turbulent parameters adopted in this paper are listed in 
Table 1 . C s , Ck,C*, C R ar e the same in the case of the solar 
model (Zhang & Li 2012a). The diffusion parameters Cti and 
C e i are set to zero, because the calculation doesn't converge 
when the diffusions of u' r T' and T'T' are taken into account. 
The turbulent dissipation parameter a = 0.85 is adopted. This 
value is due to the solar calibration based on the input physics 
and other TCM parameters. 

The age and chemical composition of each component of 
HY Vir are assumed to be the same. Interactions between the 
primary and the secondary are ignored. The stellar models 
evolve from the ZAMS to the age when the radius of the pri- 
mary stellar model fits the observation. The time step is no 
more than 0.5% of the age. All stellar models comprise more 
than 2000 mesh points. 

3. CALIBRATIONS OF STELLAR MODELS OF HY VIR 

The main aim of this paper is to test whether the solar turbu- 
lent parameters (i.e., the TCM parameters in Table 1 and the 
overshooting mixing parameter Cx = 10~ 10 ) is suitable for 
the core overshooting region and reproduce the observations 
of HY Vir. Mathematically, only one turbulent parameter can 
be derived by calibrating the effective temperatures T e / / and 
radii R of stellar models for HY Vir. I choose the most sen- 
sitive parameter, i.e., Cx, to be adjustable. The strength of 
the overshooting mixing is determined by the diffusion coef- 
ficient, which is proportional to Cx and y/k. Besides Cx, 
the TCM parameters indirectly affect the diffusion coefficient 
because of their effects on k. However, in an acceptable range 
of the TCM parameters, k doesn't change too much. 

The adjustable parameters in modeling the components of 
HY Vir are as follows: the overshooting mixing parameter 
Cx, the initial hydrogen abundance X, and the metal abun- 
dance Z. The observations make it possible to fix the ad- 
justable parameters (Cx, X, Z). The effective temperatures 
and the radii of the stars (i.e., T\ and Ri for the primary, T2 
and R 2 for the secondary) are the functions of the adjustable 
parameters (Cx, X, Z) and the age t. There are four free in- 
dependent variables and the same number of dependent vari- 
ables. Because the effective temperatures and the radii should 
be consistent with the observations, the adjustable parameters 



TABLE 2 

The observation s of HY Vir and the standard errors 
i Ilacy & FekelI20 1 II) . 



Parameter 


Primary 


Secondary 


M/Mq 


1.838 ±0.009 


1.404 ± 0.006 


l 9(Teff) 


3.836 ± 0.008 


3.816 ± 0.008 


R/Rq 


2.806 ± 0.008 


1.519 ± 0.008 


lg(L/L Q ) 


1.20 ±0.04 


0.58 ±0.04 



and the age can be solved in the sense of math. 

The scheme of solving the adjustable parameters is as fol- 
lows: (1) Use the trial parameters set (Cx, X, Z) to model 
the primary of HY Vir. (2) Obtain the effective tempera- 
ture T\ and the age t when the radius of the primary Ri 
evolves to the observational value. (3) Use the same pa- 
rameters set (Cx, X, Z) to model the secondary of HY Vir. 
The stellar model evolves from the ZAMS to the age t. Ob- 
tain the effective temperature and the radius, i.e., T2 and R2, 
of the stellar model of the secondary at age t. (4) Com- 
pare IgTi, lgT 2 and R 2 with the observations shown in Ta- 
ble 2. Use the Newtonian iteration method to revise the pa- 
rameters (Cx, X, Z), then go to step (1) until IgTi, lgT 2 , 
and R 2 are consistent with the observations within the ac- 
curacy of 0.001. Ri is calibrated to fix the age t, because 
the radius of the primary is the most sensitive variable, i.e., 
\dRi/dt\ > Max{\dR 2 Jdt\,\dlgTijdt\,\dlgT 2 /dt\). There- 
fore, to derive t by calibrating Ri is the most accurate way 
comparing with calibrating R 2 , T\, and T 2 . 

Using the above scheme, the adjustable parameters and 
their standard errors are worked out as follows: 

C x = (0.8 ±0.5) x 10~ 10 (5) 



X = 0.67 ±0.03 (6) 



Z = 0.031 ±0.007 (7) 

For the calibrated models, the age of the system is t = 
lAGyr. Those are consistent with Lacy & Fekel's(201 1) re- 
sults, which show Z = 0.027, AY/AZ = 2.0 (Y is the initial 
helium abundance) and t = (1.35 ± 0.1 )Gyr. 

The standard errors a of the parameters (Cx, X, Z) are 
calculated according to the Gaussian distribution: 

where (xi, x 2 , x a , £4, x 5 , x e ) = 

(lgT A ,lgTB,RA,RB,MA,M B ) and (2/1,2/2,2/3) 
(Cx,X, Z). The values of cr{xj) are taken from Table 
2. Although the observations of IgT, R and M are very 
accurate, Cx shows a relatively large 0. The standard errors 
of the parameters (Cx, X, Z) are mainly due to a(lgTs), 
because the contributions of j = 2 terms in Eq.(8) are 
(64%, 44%, 85%) for i= (1,2,3). 

Figure 1 shows the evolutional tracks of stellar models 
for both the primary and the secondary, with X = 0.67, 
Z = 0.031 and different Cx- Three cases are shown: Cx = 
(no overshooting mixing), Cx — 0.8 x 10 -10 (the result of 
calibration) and Cx — 1.0 x 10 -10 (suitable for downward 
overshooting region of the convective envelope in low-mass 
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3.85 |gj 3.80 




FIG. 1. — The evolutionary tracks of two components of HY Vir with 
different C'x in the diagram of (a) radius & (b) luminosity vs. effective 
temperature. Note that the temperature, luminosity and radius scales are 
different for the two components of HY Vir. Solid lines correspond to 
Cx = 0.8 X 10~ 10 , the calibrated model. Dashed lines correspond to 
Cx = 1-0 X 10 — 10 , the value for the solar case. Dotted lines correspond to 
the stellar models without overshooting mixing. The la error bars denoted 
as 'AO' and 'BO' are based on Table 2. 



stars dZhang & Lill2012rj IZhang 112012b ). It is found that the 
two cases including the overshooting mixing reproduce the 
observational effective temperatures and radii in la. Larger 
Cx results in higher effective temperature and higher lumi- 
nosity at the termination of the evolution tracks. This is be- 
cause the overshooting mixing refreshes the hydrogen content 
of the nuclear burning core and boosts the nuclear reaction in 
the core. 

4. PROPERTIES OF THE STELLAR MODELS 

4.1. The profile of chemical abundance 

The main effect of the diffusive mixing of the convec- 
tive overshooting is the modification on the profile of hydro- 
gen abundance. Figure 2 shows hydrogen abundance pro- 
files of the calibrated stellar models. Stellar models with 
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FIG. 2. — The profiles of hydrogen abundance in the stellar interior of (a) 
the primary stellar models with Xq = 0.32 and Xq = 0.50, and (b) the 
secondary stellar models with Xq = 0.55 and Xq = 0.60. The e-fold 
location of the chemical profile is defined as the location where X = Xg + 
(X(; — X$ ) /e in which Xg is the hydrogen abundance at the stellar surface. 
The solid lines correspond to the calibrated stellar models, and the dashed 
lines correspond to stellar models with 0.3Hp fully mixing in overshooting 
region for comparison. The arrows indicate the convective boundary. 



Xc = 0.50, 0.32 for the primary and stellar models with 
Xc = 0.60, 0.55 for the secondary are shown, where Xc is 
the hydrogen abundance of the convective core. The hydrogen 
abundance profiles of the stellar models with complete mixing 
in 0.3Hp overshooting region are also shown for comparison. 
The standard convective boundaries are denoted as the arrows 
in the figure. In the stellar models with complete mixing, the 
profile is discontinuous at the convective boundary when the 
convective core expands (e.g., Fig. 2(b)) or with high gradient 
when the convective core slowly contracts (e.g., Fig. 2(a)). In 
the stellar models with the diffusive mixing, X continuously 
changes from the convective boundary to the outer envelope. 

The most important property is the length of the region, 
in which chemical abundance changes due to the overshoot- 
ing mixing, since the length indicates the strength of the re- 
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freshment of hydrogen in the nuclear burning core. Because 
the profile of hydrogen abundance continuously changes, it is 
helpful to discuss by using the concept, i.e., the 'e-folding' 
length. In the overshooting region, the initial hydrogen abun- 
dance can be assumed as the surface hydrogen abundance X$, 
and the equilibrium abundance is Xc- The e-folding length 
of the chemical abundance changing region is the distance 
from the convective boundary to the e-folding location de- 
fined as X = Xs + (Xc — Xs)/e. Figure 2 shows the e- 
folding length of corresponding stellar models of two com- 
ponents. It is found that the e-folding length increases dur- 
ing the stellar evolution. In order to understand the relation 
between the e-folding length and the age of stellar model, I 
estimate the e-folding length. In the concerned evolutionary 
stage of the components of HY Vir, the stars are in main se- 
quence phase, and the stellar structures are relatively stable 
(i.e., no rapid core contraction, expansion, ignition of new el- 
ements etc.). The convective boundaries show only a little 
variations. Therefore, one can approximately study the diffu- 
sion in a 'quiescent' scene. The time required by the matter 
diffusing from the convective boundary to the e-fold location 
should be comparable with the age of the stellar model: 



t 



;2 



AD 



(9) 



ov 



where lg is the e-folding length. The diffusion time t 
l\j '(ADov) is based on the fundamental solution (i.e., the 
Green's function) of the basic diffusion equation. Defining 
Ie = ctp, Hp, and according to the asymptotical solution of 
the TCM (Zha ng & Lill2012bl) . one finds that the diffusion co- 
efficient at the e-folding location is as follow: 



D. 



ov 



-9a E /2) 



(10) 



where kc is the turbulent kinetic energy at the boundary of the 
convective core, 8 = dlnk/dlnP is the exponential decreas- 
ing index of the turbulent kinetic energy in the overshooting 
region. According to the asymptotical solution, the parame- 
ters of the TCM adopted in this paper leads to 9 = 4.8. There- 
fore, Eq.(9) can be rewritten as: 



2 



12.6 



Cx,~wWcitG 



H 



P.10 



(11) 



where C x ._ 10 = Cx/ICT 10 , W CA = ^/(lO^crus- 1 ), 
t G = t/(lGy) and H P , W = H P /\lO w cm). Equation (11) 
indicates that the e-folding length increases during the stellar 
evolution, which has been found in Fig. 2. 

Figure 3 shows the comparison between the e-folding 
length predicted by Eq.(l 1) and the numerical results of stel- 
lar models. It is found that the predicted e-folding length is 
consistent with the numerical results. The predicted e-folding 
length is zero when t — > 0. However, the e-folding length in 
stellar models is not zero when t — > 0. This is because the 
nuclear burning outside the convective core results in a chem- 
ical abundance changing region. Equation (11) takes only 
the overshooting mixing into account, and ignores the nuclear 
burning. It should be emphasized again that Eq.(ll) holds 
only if the stellar structure is relatively stable. It is invalid 
for stellar models after the end of the main sequence phase, 
during which the convective core starts to contract rapidly. 

4.2. The temperature gradient in the overshooting region 
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FIG. 3. — The e-fold length parameter ole = Ie/Hp of the chemical 
abundance changing region during the stellar evolution. The solid lines are 
corresponding to the stellar model, and the dashed lines are corresponding to 
the prediction by Eq.(l 1). 



In this paper, the stellar structure equations are solved to- 
gether with the TCM. In the convection zone, the MLT, which 
is used in the standard stellar evolutionary calculations, is re- 
placed by the TCM. Comparing with the MLT, the localized 
TCM shows similar results of temperature gradient in the con- 
vection z one. H owever, the TCM results in non-zero turbulent 
heat flux u' r T' in the overshooting region. The temperature 
gradient in the overshooting region is modified according to 
Eq.(4). 

Figure 4 shows the temperature gradient near the convective 
boundaries. It can be found in Fig. 4 that V > W p in the over- 
shooting region. This is a c ommon property of non-local tur- 
bulent convection the ories (|Xiong|[l985[ 1 1 989t ICanutoll 199 
iXiong & Deng l200ll iDeng et al.1 120061: iZhang & Lil 1200 " 
These theories show u' r T' < in the overshooting region be- 
cause the buoyancy prevents turbulent motions, thus V > V.r 
according to Eq.(4). 

It is found in Fig.4 that, in the case of high Peclet number at 
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b) The secondary 



b) The primary, X =0.32 




0.44 



0.42 



0.40 



0.38 



0.36 



39.5 



FIG. 4. — The temperature gradient near the convective boundaries. The 
stellar models are (a) the primary stellar model with Xq = 0.32 and (b) the 
secondary stellar model with Xq = 0.55. The left part of each sub-figure 
is for the base of the convective envelope, and the right part is for the bound- 
ary of the convective core. P e (C) is the Peclet number at the convective 
boundary. l at { is the length of the adiabatic overshooting region. 

the convective boundary, i.e., P e (C) 3> 1, there is a small adi- 
abatic overshooting region adj oining the convecti ve boundary. 
This has been pointed out by Zhang & Li (2012b) in their the- 
oretical analysis that, in the overshooting region with P e 3> 1, 
Cei = (adopted in this paper) results in an adiabatic over- 
shooting region. The length of the adiabatic overshooting 
region l a d can be e stimated by using the parameters of the 
TCM (Zhang & Lill2012bT) . The parameters adopted in this 
paper shows l a d ~ O.OlHp. The diffusion of the turbulent 
temperature fluctuation V = T'T' is absent because of some 
problems of numerical calculations. If it is taken into account, 
there should be no adiabatic overshooting region. In most part 
of the overshooting region even with P e {C) ^> 1, the temper- 
ature gradient is close to the radiative one. This may be for the 
reason that, in the overshooting region with P e (C) 3> 1, the 
correlativity o f turbulent velocity and temperature is almost 
equal to zero (Xiona [l985l : IXi ong & Deng 2001: lDeng etalj 




16 |nP32 



FIG. 5. — The turbulent variables of the stellar models of the primary, k is 
turbulent kinetic energy, U 



T'T' the 



u' r T' the turbulent heat flux, V 
turbulent temperature fluctuation, P e the Peclet number, uj = u' r u' r /(2k) 
the anisotropic degree. The arrows indicate the convective boundary. '+' and 
'-' around the arrows indicate the convection zone and overshooting region, 
respectively. 



I2006t IZhang &XHl2012bh . thus turbulent motions can hardly 
transport heat. 

Although the TCM modifies the temperature gradient in the 
overshooting region, the modification is small. Zhang & Li's 
(2012a) results have indicated that, comparing with the over- 
shooting mixing, the temperature gradient modification based 
on the adopted TCM parameters is ignorable on affecting the 
stellar structure and evolution. 

4.3. The turbulent variables in the overshooting region 
The turbulent variables involved in the TCM are as follows: 



u' r u' r is the radial turbulent kinetic energy, k total turbulent 
kinetic energy, u' r T' the turbulent heat flux, T'T' the turbulent 
temperature fluctuation. 

The overshooting region can be classified into two kinds: 
the high P e overshooting region with P e 3> 1 and the low 
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The secondary, X c =0.55 




i i . . . 



10 15 20 25 | n p 30 35 40 



TABLE 3 

The properties of the turbulent variables in the high P e 

OVERSHOOTING REGION: NUMERICAL RESULTS VS . THEORETICAL 
PREDICTIONS. 





C Fig.5(a) 


C Fig.5(b) 


C Fig.6 


E Fig.6 


Asy. 


Vkc 


2060 


2760 


1410 


16500 




v^CTMOD) 


2120 


2510 


1630 


15200 




dlnk/dlnP 


4.5 


4.4 


4.6 


-4.8 


±4.8 


dln(-u' r T')/dlnP 


6.8 


7.1 


6.7 


-7.2 


±7.2 


dln(T 7 T 7 )/dlnP 


4.6 


4.5 


4.5 


-4.6 


±4.8 


u = «/(2fc) 


0.28 


0.28 


0.28 


0.28 


0.28 



NOTE. — kc is the turbulent kinetic energy at the convective bound- 
ary. The 'Asy.' column show the results predicted by the asymptotical 
analysisi Zhang & Li 20123). 'V^C (TMOD)' is the value of y/kg p redicted 
by the method called 'the maximum of diffusion' (Zhang & Li 2012b). The 'C 
Fig. 5(a)' column means the core overshooting in Fig. 5(a), and the 'E Fig.6' 
column means the envelope overshooting in Fig.6. 



FIG. 6. — Similar to Fig. 5, but for the stellar models of the secondary. 



P e one. The core overshooting region is always the high P e 
one because of the short free path of photons in the stellar in- 
terior. The type of the downward overshooting region of the 
convective envelope depends on the temperature at the con- 
vective boundary Tbce- Commonly, it is the high P e over- 
shooting region for IgTscE > 5.5 or the low P e one for 
IqTbce < 5.5. The properties of the two kinds of over- 
shooti ng region have been found to be different with each 



other (X iong & Chenll992tlXiong & DengHOOltlZhang &H 



l2012bl) . The properties of the high P e overshooting region in 
th e framework of Li & Yang's (2007) TCM have been studied 
bv lZha ng & Li ( 2012bl) . However, the property of the low P e 
overshooting region is not as clear as the high P e one. 

Figures 5-6 show the turbulent variables in the stellar inte- 
rior. For the primary, the ZAMS model and the model with 
Xc = 0.32 are shown. For the secondary, only the stellar 
model with Xc — 0.55 is shown, because the structure of the 
turbulent variables shows only a little variations during the 
concerned evolutionary range. At all convective boundaries, 
u' r T' and T'T 1 are equal to zero because their diffusions are 
ignored in the calculations. 

As shown in Figs. 5-6, T'T' in the core overshooting region 
is larg er than in the core convection zone. This is coincident 
with IXiongl (119851) . The turbulent temperature fluctuation 
T'T' represents the difference of the temperature between the 
turbulent flows and the environment, i.e., AT. There is al- 
ways P e ^> 1 in the core convection zone. In the P e ^> 1 
convection zone, the turbulent heat transport is very efficient, 
leading to V w V ' ad . Therefore, the temperature of the turbu- 
lent flows is almost equal to the temperature of the environ- 
ment, i.e., AT w 0, thus T'T' w in the P e > 1 convection 
zone. In the overshooting region, however, V is close to V n 
so that AT is significant. 

In the high P e overshooting region (in the right parts of 
Fig.5(a-b), and in Fig.6), the turbulent variables k, u' r u' r , 
T'T' decrease exponentially, and the anisotropic degree lu = 
u' r u' r /(2k) almost doesn't change. Thos e are consistent with 
the asymptotical analysis of the TCM dZhang & Lill2012bl) 
and Xiong 's non-local turbulent model (see, e.g. , [Xiong 
( fl98ll[T989l) : IXiong & Dend d200lb : lDeng & Xiond d2008l) )7 



The overshooting region between 
two convection zones 

k 05 /10 5 

U/10 6 

V/2000 / 

10 4 F/F 




11.0 



11.5 



InP 



12.0 



12.5 



FIG. 7. — The turbulent variables in the overshooting region which lays 
between two convection zones. The dotted line is the ratio of the convective 
flux (see text) to the total flux Fq/F = 1 — V/Vjj. The stellar model is 
the ZAMS of the primary. Symbols are the same as Fig. 5. 



Table 3 compares the properties of the turbulent variables in 
the high P e ove rshooting region wi th the results of the theo- 
retical analysis (Zhang & Li 2012b). It is found that the nu- 
merical results are in agreement with the theoretical analysis. 
The high P e overshooting region extends to the location of 
P e ~ 1 and then becomes the low P e one. 

In the low P e overshooting region, the thermal dissipation 
dominates, thus the turbulent variables k, u' r u' r , T'T' decrease 
super-exponentially and u decreas es to zero accordin g to the 
properties of the TCM equations (Zhang & Li 2012b). Those 
properties can be found in Fig. 5. The overshooting is finally 
cut-off in the region P e -C 1. According to the definition of 
P e , P e oc y/k, thus P e exponentially decreases in the over- 
shooting region with the index 9/2. As a consequence, the 
distance (in Hp) from the convective boundary to the cut-off 
location depends on P e (C). Larger P e (C) leads to longer dis- 
tance to the cut-off location. This can be found in Figs. (4-6) 
by comparing P e (C) and the the distance from the convective 
boundary to the cut-off location of each overshooting region. 
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There is a thin overshooting region laying between two con- 
vective shells, which are caused by H and Hell ionization re- 
spectively, in the ZAMS stellar model of the primary as shown 
in Fig. 5a. The convective shells merge together during the 
stellar evolution, thus this overshooting region disappears in 
the Xc = 0.32 stellar model. The details of the turbulent 
variables in the thin overshooting region are shown in Fig. 7. 
The turbulent velocity \fk ~ 10 5 cm/ s in the region is larger 
than in the convective core, so that two convective envelopes 
can be considered to be dynamically connected. However, the 
minimum of the ratio of the convective flux to the total flux is 
about \F C /F\ = |l-V/V fl | ~ 10" 5 , where F c = pc P u' r T' 
is the convective flux. Moreover, P e < 1 in the two con- 
vective shells and the overshooting region. These conditions 
prevent the turbulent heat transport between two convective 
shells. The convective envelopes can therefore be considered 
to be thermally disconnected. 

5. CONCLUSION AND DISCUSSION 

In this work, the turbulent convection model (TCM) is in- 
volved in the stellar structure and evolution in order to study 
the convective core overshooting in the stellar models of 
eclipsing binary star HY Vir. Using new accurate observa- 
tions, I calibrate the stellar effective temperatures and radii 
of two components of HY Vir, and then obtain the required 
overshooting mixing parameter Cx = (0.8 ± 0.5) x 10 -10 . 
The turbulent parameters (i.e., the TCM parameters and the 
overshooting mixing parameter Cx), which is proper for the 
downward overshooting region of the convective envelope in 
low-mass stars, can reproduce the observational radii and ef- 
fective temperatures of the two components of HY Vir. It 
indicates that the diffusive overshooting mixing based on the 
turbulent velocity described by the TCM can also be applied 
to the core overshooting. This result encourages us to apply 
the TCM to the convective overshooting. 

The diffusive overshooting mixing causes a continuous pro- 
file of the hydrogen abundance. The e-folding length of the 
region, in which the chemical abundance changes due to the 
overshooting mixing, can be estimated by using Eq.(ll). It 
is found that the e-folding length increases during the stellar 
evolution. 

The overshooting region can be classified into two kinds: 
the high P e overshooting region in which P e 1 and the 
low P e one. The properties of the high P e oversho oting region 
are in agreement with the theoretical analysis (Zha ng & Lil 
l2012bl) . The turbulent variables (turbulent kinetic energy k, 
turbulent heat flux u' r T' and turbulent temperature fluctua- 
tion T'T') decrease exponentially in the high P e overshoot- 
ing region. The properties of the low P e overshooting re- 
gion are very different from the high P e one. In the low P e 



overshooting region, the turbulent variables decrease super- 
exponentially and are finally cut-off. 

The diffusions of u' r T' and T'T' are not taken into ac- 
count in this paper because of some problems of numerical 
calculations. However, Cx of the calibrated stellar models 
for HY Vir should not change significantly if those diffusions 
are present. Those diffusions may affect the calibrated results 
mainly on two ways: (i) modifying the turbulent heat flux and 
(ii) affecting the core overshooting mixing which plays an im- 
portant role on the stellar evolution. Although those diffusions 
can redistribute u' r T' in the stellar interior, they should't sig- 
nificantly change the integral value of u' r T'. The effective 
temperature of stellar mod els sh ould be insensitive to those 
diffusio ns. The diffusion of u' r T' is ignorable in the region of 
P e 3> 1 (Zhang & Li 2012b), thus it doesn't affect the proper- 
ties of the core convection and overshooting. The diffusion of 
T'T' can affects the exponential index of the turbulent kinetic 
energy 6(— dlnk/dlnP). The solar value for the diffusion of 
TT is small that C el = 0.02 dZhang & Lill2012ab . It should 
modify the turbulent properties in the core overshooting re- 
gion slightly. Cx of the calibrated model should change a 
little if the diffusion of T'T' is present. However, since the 
relative standard error of Cx is about 60%, the solar value 
Cx = 10~ lf} should also reproduce the required radius and 
effective temperature of HY Vir in la. 

In other treatments of overshooting with incomplete mix- 
ing mentioned in Section 2.1, there is no information on 
how the temperature gradient profile is modified near the 
convective boundary. The helioseismic investigation by 
IChristensen-Dalsgaard et aD (1201 ll) can't judge which of 
them is suitable for the solar case. However, the mixing can 
affect the sound speed profile. It is well known that, in the 
standard solar model, there is a bump of sound speed dif- 
ference below the solar convection zone. In order to elim- 
inate the bump, the proper diffusion coefficient at the base 
of the solar convection zone is on the magnitude order of 
10f_~ 10 3 (see, e.g.. IChristensen-Da lsgaard & Di Mauro 
d2007l) : IZhang & Lil d2012al) ). The mixing should be effec- 
tive in about O.li?, which is the width of the bump, below the 
base of convection zone. The helioseismic inversion requires 
the diffusion coefficient to satisfy those conditions in the solar 
case. 
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